Hole-doped, High- Temperature Superconductors Li x BC, Na x BC and C x : A 

Coherent-Potential-Based Prediction 
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Using density-functional-based methods, we show that the hole-doped Li x BC and Na x BC in 
PQz/mmc crystal structure and C x in graphite structure are capable of showing superconductivity, 
possibly with a T c much higher than that of MgB-2. We use full-potential methods to obtain the 
optimized lattice constants a and c, coherent-potential approximation to describe the effects of disor- 
der, Gaspari-Gyorffy formalism to obtain the electron-phonon coupling constant A, and Allen-Dynes 
equation to calculate T c as a function of hole concentration in these alloys. 

PACS numbers: 74.25.Jb, 74.70.Ad 



The search for new MgB2-like superconductors |jj has 
proved to be elusive so far, in spite of intense experimen- 
tal effort. The difficulty in finding new materials with 
superconducting properties similar to that of MgB 2 @-||] 
may partly be due to the delicate nature of interaction 
responsible for superconductivity (SC) In MgB 2 , 

the relatively high superconducting transition tempera- 
ture, T c , of 39 K is achieved through a coming together 
of various aspects of phonon-mediated electron-electron 
interaction, such as (i) strong in-plane B — B covalent 
bonding , (h) electron and hole-like cylindrical Fermi 
sheets along r to A > (iii) a strong coupling of the a 

states to the in-plane B — B bond-stretching mode |^,||j8| 



including an-harmonic effects 17U8|], (iv) B p x ^ contribu- 
tion to the total density of states (DOS) at the Fermi 
energy, Ep, and the total DOS at Ep [B](J, and others. 

Most of the experimental as well as theoretical efforts 
have used layered and metallic materials as the start- 
ing point to obtain materials with above-mentioned elec- 
tronic properties leading to high-temperature SC but 
without success. It may be that we have to follow a dif- 
ferent route to obtain new MgB 2 -\ike superconductors 
with high T c . Such an alternative route is provided by 
hole-doping of the layered but semiconducting or insulat- 
ing materials (9). 

With the aim of developing materials with MgBz-iike 
electronic properties, we have studied how hole-doping 
can be used to tailor the electronic structure of Li x BC, 
Na x BC and C x alloys. In particular, it turns out that 
the hole-doping makes Li x BC, Na x BC and C x alloys 
similar to MgB 2 in more ways than one, making these al- 
loys promising candidates for showing high-temperature 
SC. 

In the following, we present the results of (i) fully- 
relaxed, full-potential electronic structure calculations 
for LiBC and NaBC in P6 3 /mmc crystal structure and 
C in graphite structure, using density-functional-based 
methods, (ii) charge- self-consistent, coherent-potential- 
based electronic structure calculations for the hole-doped 



Li x BC, Na x BC and C x alloys, and (iii) the calculations 
of T c in hole-doped Li x BC, Na x BC and C x alloys with 
x ranging from 1 to 0.2 (hole concentration is equal to 
(1-x) ). 

For full-potential calculations, we have used ABINIT 
code @, based on psuedopotentials and planewaves, 
and optimized the lattice constants a and c of LiBC, 
NaBC and C. The effects of hole doping in Li x BC, 
Na x BC and C x alloys have been studied using Korringa- 
Kohn-Rostoker coherent-potential approximation in the 
atomic-sphere approximation (KKR-ASA CPA) method 
]TT] , p~2[ . For calculating T c as a function of x in Li x BC, 
Na x BC and C x alloys, we have used the Gaspari-Gyorffy 
formalism |D| to evaluate the electron-phonon (EP) cou- 
pling constant A, and Allen-Dynes equation H] to cal- 
culate T c . 

Based on our calculations, described below, we find 
that the hole-doped Li x BC and Na x BC and C x are 
capable of showing SC, possibly with high as 

- 105 if for Li . 5 BC ||, - 106 if for Na . 6 BC and 
~ 60 K for Co. 4. Such a high T c in these alloys be- 
comes possible partly because the hole-doping modifies 
the electronic structure, discussed in terms of the density 
of states and spectral function along r to A direction, and 
brings it closer to that of MgB 2 . 

Before describing our results in detail, we provide some 
of the computational details. The structural relaxation 
of LiBC, NaBC and C was carried out by the molecu- 
lar dynamics program ABINIT as described in Ref. ||. 
The charge self-consistent electronic structure of Li x BC, 
Na x BC and C x alloys as a function of x has been calcu- 
lated using the KKR-ASA CPA method. We have used 
the CPA rather than a rigid-band model because CPA 
has been found to reliably describe the effects of disor- 
der in metallic alloys 0,Q, including MgB 2 alloys [ fl5| |. 
We parametrized the exchange-correlation potential as 
suggested by Perdew-Wang |l6| within the generalized 
gradient approximation. The Brillouin zone (BZ) inte- 
grations for charge self-consistency and the DOS calcu- 
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FIG. 1. The total density of states calculated at the opti- 
mized lattice constants using the ABINIT program. For com- 
parison the DOS for MgBi has been doubled. The vertical 
dashed line indicates Ef- 



TABLE I. The calculated lattice constants a and c in 
atomic units. The experimental lattice constants are shown 
in the parentheses. 



Compound 


a 


c 


LiBC 


5.16 (5.202) 


13.24 (13.342) 


NaBC 


5.25 


15.62 


C 


4.63 (4.650) 


12.62 (12.649) 


MgB 2 


5.76 (5.834) 


6.59 (6.657) 



lations were carried out using 1215 and 4900 k-points in 
the irreducible part of the BZ, respectively. For DOS 
(spectral function) calculation, we added a small imagi- 
nary component of 1 (2) mRy to the energy. The lattice 
constants for Li x BC, Na x BC and C x alloys as a func- 
tion of x were kept fixed to the values given in Table I. 
To reduce ASA-related errors in open structures we in- 
troduced empty spheres. The Wigner- Seitz radii for Li 
and Na were slightly larger than that of B and C. The 
maximum I value used was l m ax = 2. 

The A was calculated using Gaspari-Gyorffy formal- 
ism with the potentials of Li x BC, Na x BC and C x ob- 
tained with the KKR-ASA CPA method. Subsequently, 
the variation of T c with x was calculated using Allen- 
Dynes equation. The average value of phonon frequency 
coin for MgB-i was taken from Ref. ^ , while the value of 
fi* was kept fixed at 0.09. 

First, we discuss the changes brought about by hole- 
doping in the various aspects of the electronic structure of 
Li x BC, Na x BC and C x alloys in terms of (i) total DOS, 
(ii) partial B and C p-contributions to the total DOS at 
Ef, and (iii) spectral function along T to A evaluated at 
Ef- Then we describe the calculated T c for these alloys 
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FIG. 2. The total DOS (solid line) for Li x BC (upper 
panel), Na x BC (middle panel), and C (lower panel) at 
x = 0.5, calculated using the KKR-ASA CPA method. For 
comparison the total DOS for MgB^ (dashed line) is also 
shown. The vertical dashed line indicates Ef- 

as a function of hole concentration, which is the main 
result of the present work. 

We show in Fig. 1 the total DOS for LiBC, NaBC 
and C alloys calculated using the ABINIT program at 
the optimized lattice constants as given in Table I. For 
comparison we also show the DOS for MgB2 obtained as 
described in Ref. [|| . The compound NaBC has a higher 
cj a than LiBC primarily because of the larger Na atoms. 
The difference in c/a ratio leads to significant changes 
between the electronic structure of LiBC and NaBC. 
By comparing the DOS of LiBC, NaBC and C with 
MgB-2, it is clear that if Ep can be moved suitably inside 
the valence band then the DOS at Ep in these alloys can 
be made comparable to the DOS at Ef in MgB%. One 
of the ways of moving the Ef inside the valence band is 
through hole-doping, as described next. 

The changes in the electronic structure of Li x BC, 
Na x BC and C x alloys as a function of x has been studied 
for x ranging from 1 to 0.2. We find that with increasing 
hole concentration Ep moves inside, as expected, and 
consequently these alloys become metallic. As an ex- 
ample, we show in Fig. 2 the total DOS of Li^^BC, 
Na , 5 BC and Co. 5 alloys, calculated as described earlier. 
From a comparison with the DOS of MgBi, calculated 
using the same approach and also shown in Fig. 2, it is 
clear that the DOS around Ef in Li^.^BC, Nao.^BC and 
Co. 5 alloys are comparable to that of MgB2- Our calcula- 
tion shows that the hole concentrations of 50% — 80% are 
ideal for matching the total DOS of Li x BC , Na x BC and 
C x alloys with that of MgB 2 - Knowing the importance 
of the total DOS as well as the in-plane Bp-contribution 
to the total DOS at Ep in MgB 2 , the alloys Li x BC, 
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FIG. 3. The p x t v ) component of B (solid line) and C 
(dashed line) partial DOS around Ef for (a) Lio.sBC, (b) 
Na . 5 BC, and (c) C . 5, calculated using the KKR-ASA CPA 
method. For comparison the p x ( y ) component of B partial 
DOS in MgBi (dotted line) is also shown. The vertical dashed 
line indicates Ef- 



FIG. 4. The calculated spectral function along F to A di- 
rection , evaluated at Ef, for x — 0.3 (solid line) and x = 0.5 
(dashed line) in (a) Li x BC, (b) Na x BC, and (c) C x alloys. 
For comparison, the spectral function for MgB2 (dotted line) 
is also shown. In (b) the spectral function for Nao.sBC is out 
of range of the plot. 



Na x BC and C x will be good candidates for showing SC 
only if the B and Cp-contributions to the total DOS at 
Ep are comparable to that of MgB%. 

In Fig. 3 we compare the p x {y) contributions of B 
and C to the total DOS around Ep with the correspond- 
ing contribution from MgB-2- The Bp x r v \ contribution 
in Liq.^BC is almost identical to that of MgB2, while 
the C p x (y) contribution is slightly higher, indicating a 
strong B — C covalent bonding. In Na^.^BC the Bp x t y \ 
contribution is slightly lower than that in MgB 2 . Thus 
a further increase in hole concentration will enhance the 
p x ( y ) contribution to the total DOS. On the other hand 
in the case of Co. 5 the C p x ( v ) DOS is identical to that of 
the Bp x ( y } contribution in MgB 2 . 

For Li x BC, Na x BC and C x alloys to remain in the 
contention for showing M gB 2 -\ike superconducting prop- 
erties, a favorable matching of the total DOS and the 
partial Bp x i y \ and/or C p x (y) contributions to the total 
DOS at Ef must be followed by similarities in the Fermi 
surfaces of these hole-doped alloys with that of MgB 2 , 
especially along r to A. A measure of this similarity is 
provided by comparing the spectral function, A(k, Ep), 
of MgB 2 with Li x BC, Na x BC and C x alloys along T to 
A [jl5|| . Such a comparison is shown in Fig. 4, where we 
have plotted A(k, E F ) along T to A for Li x BC, Na x BC 
and C x alloys with x = 0.3 and 0.5, as well as for MgB 2 . 
We find that A(k, Ep) for Lio^BC is closest to that of 
MgB 2 , indicating similarities in the Fermi surfaces of the 
two compounds along T to A. In C x the spectral func- 
tions for x = 0.3 and 0.5 are 7 — 8 times larger than 
the corresponding values in MgB 2 , while in Na^ ^BC 



it is 40 — 45 times larger. Based on the comparison of 
the various aspects of the electronic structure of Li x BC, 
Na x BC and C x alloys with MgB 2 , we find Li x BC to be 
the most promising candidate for showing MgB 2 -like SC 
followed by C x . Such a conclusion is borne out by our ex- 
plicit calculation of T c as a function of hole concentration 
in Li x BC, Na x BC and C x alloys. 

The calculated T c as a function of hole concentration in 
Li x BC, Na x BC and C x alloys is shown in Fig. 5, where 
we also show the variation in T c in MgB 2 . Clearly, hole- 
doped Li x BC, expected to have a T c of around 105 K 
with 50% hole concentration, offers the best hope for im- 
proving the T c of MgB 2 -type superconductors. The EP 
coupling is much stronger in Li x BC (A ~ 1.7 for x = 0.5) 
than in MgB 2 or any other alloys that we have consid- 
ered. Our results for Li x BC are consistent with the re- 
sults of Ref. ||. Our calculation also shows that Na x BC, 
if it can be synthesized, will also have a T c comparable to 
that of Li x BC albeit at a higher hole concentration. The 
change in T c of MgB 2 with increasing hole concentration 
is not large. Surprisingly, C x is also found to show SC at 
elevated temperature. 

The work presented here can be further improve d by 
optimizing the lattice constants a and possibly c [p.7[ , 
within the CPA, as a function of hole concentration. 
However, the results presented in this paper are very ro- 
bust and, in our opinion, a slight change in lattice con- 
stants due to optimization will not change the main con- 
clusions of the present work. 

In conclusion, we have presented the results of a fully- 
relaxed, full-potential electronic structure calculations 
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FIG. 5. The calculated variation in T c as a function of hole 
concentration for Li x BC (filled square), Na x BC (filled tri- 
angle), C x (star), and Mg x B2 (filled circle). 
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for LiBC and NaBC in PQ^/mmc crystal structure and 
C in graphite structure. Using CPA, we have discussed 
the effects of hole-doping in Li x BC, Na x BC and C x al- 
loys in terms of (i) total DOS as well as Bp contribution 
to it and, (ii) the spectral function along T to A, and 
shown that the electronic structure of the hole-doped al- 
loys is similar to that of MgB^. Our calculated T c clearly 
shows Li x BC and C x to be superconducting in the range 
0.2 < x < 0.8, and Na x BC to be superconducting in the 
range 0.2 < x < 0.65, with Li .$BC showing a possible 
T c of 105 K. 
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